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Abstract 

We describe recent extensions of the program SPheno including flavour aspects, CP-phases, R-parity 
violation and low energy observables. In case of flavour mixing all masses of supersymmetric particles are 
calculated including the complete flavour structure and all possible CP-phases at the 1-loop level. We 
give details on implemented seesaw models, low energy observables and the corresponding extension of 
the SUSY Les Houches Accord. Moreover, we comment on the possiblities to include MSSM extensions 
in SPheno. 

1. Program Summary 
Program title: SPheno 

Program Obtainable from: http://projects.hepforge.org/spheno/ 
Programming language: F95 

Computers for which the program has been designed: PC running under Linux, should run in 
every Unix environment 
Operating systems: Linux, Unix 

Keywords: Supersymmetry, renormalization group equations, mass spectra of supersymmetric models, 
Runge-Kutta, decays of supersymmetric particles, production 

Nature of problem: The first issue is the determination of the masses and couplings of supersymmetric 
particles in various supersymmetric models the R-parity conserved MSSM with generation mixing and 
including CP-violating phases, various seesaw extensions of the MSSM and the MSSM with bilinear R- 
parity breaking. Low energy data on Standard Model fermion masses, gauge couplings and electroweak 
gauge boson masses serve as constraints. Radiative corrections from supersymmetric particles to these 
inputs must be calculated. Theoretical constraints on the soft SUSY breaking parameters from a high 
scale theory are imposed and the parameters at the electroweak scale are obtained from the high scale 
parameters by evaluating the corresponding renormalization group equations. These parameters must 
be consistent with the requirement of correct electroweak symmetry breaking. The second issue is to use 
the obtained masses and couplings for calculating decay widths and branching ratios of supersymmetric 
particles as well as the cross sections for theses particles in electron positron annihilation. The third issue 
is to calculate low energy constraints in the B-meson sector such as BR(6 sj), AAIb^, rare lepton 
decays, such as BR(/x — > 67), the SUSY contributions to anomalous magnetic moments and electric dipole 
moments of leptons, the SUSY contributions to the rho parameter as well as lepton flavour violating Z 
decays. 

Solution method: The renormalization connecting a high scale and the electroweak scale is calculated 

by the Runge-Kutta method. Iteration provides a solution consistent with with the multi-boundary 

conditions. In case of three-body decays and for the calculation of initial state radiation Gaussian 

quadrature is used for the numerical solution of the integrals. 

Restrictions: In case of R-parity violation the cross sections are not calculated. 

Running time: 0.2 second on a Intel(R) Core(TM)2 Duo CPU T9900 with 3.06GHz 
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2. Introduction 

In its orginal version SPheno had been designed to calculate the spectrum in the MSSM neglecting 
any effects due to generation mixing and/or CP violation Moreover, the two- and three-body decays 
of the SUSY particles as well as of the Higgs bosons can be calculated as well as the production rates of 
these particles in e+e^ annihilation. The code itself is written in FORTRAN 95. 

The program has been extended to include flavour aspects, CP-violation and R-parity violation. 
Moreover, different variants of the seesaw mechanism have been implemented. In this paper we describe 
the corresponding changes and implementations. Details on the algorithms used can be found in the 
orginal manual Moreover, we give in the appendices the default values for various flags as well as the 
error coding. 

3. Extensions with MSSM particle content at the electroweak scale and conserved R-parity 

SPheno has been extended to include flavour and CP violating phases using the SLHA2 conventions 
Q for the general MSSM. For this purpose the complete flavour structures including CP-phases have 
been implemented in the RGBs at the 2-loop level using the formulas given in Q . We have extended the 
formulas of 0] for the 1-loop masses to account for the flavour structures, e.g. we calculate the 1-loop 
corrected 6x6 mass matrices for squarks and charged leptons and the 1-loop corrected 3x3 mass matrix 
for sneutrinos Moreover, we take into account all possible phases in the calculation of all mass 

matrices at the 1-loop level but for one exception: we do not consider the loop induced mixing between 
the scalar and pseudoscalar Higgs bosons. 

Beside the calculation of the spectrum also the decay routines have been extended to calculate all 
possible two- and three body decays of supersymmetric particles and Higgs bosons at leading order 
including flavour effects. In addition the set of low energy observables has been extended as described in 
section [5] 

Within the MSSM several model classes are implemented 

• High scale models like mSUGRA, GMSB, AMSB 

• A SUGRA scenario where all soft SUSY breaking parameters are given freely at the GUT scale 
which is determined usually via the condition gi(MGUT) = 52(-M'gut)- However, this scale can be 
set to a fixed value using entry 31 in block SPheno Input, see section [B. 2. 51 

• All MSSM parameters specified at the electroweak scale Qewsb with a user specified value for 
Qewsb- 

In all cases it is assumed that the required input is given via the SLHA convention 0, . 

In addition several classes of neutrino mass models have been included with additional states at high 
energy scales: 

• A seesaw I model with different masses for the right-handed neutrinos. The corresponding particle 
content can be chosen setting the entry 3 of the block MODSEL as described in section 16.1.21 The 
parameters are set using the blocks MNURIN and YNUOIN, see sections 16.2.21 and I6.2.7[ respectively. 
Using this the results of 0, Q have been obtained. 

• Two variants of seesaw II model have been implemented. Here one can either choose a pair of 
SU{2) triplets or a pair of SU{5) 15-plets to generate neutrino masses. The first version uses the 
formulas [10| including the corrections presented in JJ| and 2-loop RGBs for the gauge couplings 
and gaugino mass parameters as used in [l^]. Here the blocks Higgs3IN and YTOIN have to be used 
to transfer the data, see sections 16.2.11 and 16.2.81 This variant is faster if one uses 2-loop RGBs but 
less accurate in particular for low seesaw scales [l^ • 

In case of a pair of SU{5) 15-plets a second variant has been implemented using the complete 2-loop 
RGBs and corresponding threshold corrections at the seesaw scale as described in 14]. Here the 
blocks Higgs3IN and YTIN have to be used to transfer the data, see sections 16.2.11 and 16.2.91 
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A seesaw III model with three SU{5) matter 24-plets usmg the com plet e 2-loop RGEs and cor- 
responding threshold corrections at the seesaw scales as described in [l3|. The blocks MWMIN and 
YB3IN can be used to set the parameters, see sections [6.2.31 and 16.2.61 



• A minimal SU{5) model as described in (15[. The corresponding particle content can be chosen 
setting the entry 3 of the block MODSEL as described in section 16.1.21 The parameters are set using 
the blocks MNURIN and YNUOIN, see sections and KT7\ respectively. The additional SU{5) 

parameters can be set extending the block MINPAR as described in section lB.l.ll 

Note, that in these models the particle content at the electroweak scale is the same as in the usual MSSM 
and that the differences are only due to the modified evaluation of the parameters. 



4. R-parity violation 

Currently the bilinear model is implemented, i.e. extending the superpotential by the terms 

= e,L,Hu (1) 

and the corresponding soft SUSY breaking terms. In this class of models neutrino physics can be ex- 
plaining due to the mixing of neutralinos with neutrinos and by loop contributions. The corresponding 



details can be found in [ig, |17[ • The same parameters giving rise to neutrino masses also lead to the 
decajfjaf the LSP and, thus, there are correlations between the LSP decay properties and neutrino physics 
li[ii0,|2l|,[2i,|23|. 



Here one has two options 

• use within SLIIA2 the blocks EXTPAR, RVSNVEVIN and RVKAPPAIN to specify the model parameters 
at the electroweak scale 

• use one of the high scale models mSUGRA, GMSB or AMSB to calculate the R-parity conserving 
parameters at the electroweak scale. The R-parity parameters are then added at this scale using 
one of the two possiblities 

- add them using the blocks RVSNVEVIN and RVKAPPAIN 

— use the flag 91 of the block SPhenoInput as described in section IB. 2. 51 to calculate the and 
the sneutrino vacuum expectations values Vi such, that neutrino physics is respected. The 
corresponding neutrino data can be specified in block NeutrinoBoundsIn, see section [5.2.41 

In this class of models the mass matrices are calculated at tree- level except for the neutrino/neutralino 
mass matrix which requires the inclusion of the full 1-loop contributions. Moreover, all possible R-parity 
violating decay modes are calculated. 



5. Low energy observables 

In this section we summarize the main references from which the formulas for the corresponding 
implementation have been taken. Moreover we give implementation specific details whenever necessary. 
For the moment the low energy observables are only calculated if there is an effective MSSM at the 
electroweak scale and conserved R-parity. 

5.1. B-physics observables 

The following observables are calculated in SPheno: BR{b — > 57), BR{b — ?> sji^fi^), BR{b — > 
sT,^^i^i)^ BRiB^ M+M"), BR{B° BR{Bu t+v), AMbo and AM^o. For the cal- 

culation of the Wilson coefhcients we use running couplings and SUSY masses which are in general 
evolved at the scale Q = mz. The only exception is BR{b 37) as we use here the formula of ref. 
where the corresponding coefficients have to be given at the scale Q = 160 GeV. For the calculation of 
the Wilson coefficients and the corresponding observables we have used 
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Table 1: Parameters used in the calculation of the B-physics observables. 



T^o — 1.525 ps TgQ — 1.472 ps = 1.638 ps = 0.45 ps 
/i3 = 190MeV BBa = 1.22 /sj = 230 MeV Bs, = 1.26 
Mbo = 5.2795 GeV M^q = 5.3663 GeV ??b = 0.55 



• BR{b ^ 57) [aaliaHl; the value given is for > 1.6 GeV and rric/mb = 0.23. 

• BR{h ^ s^+^jL-) 

• BR{h ^ sY^^v.v,) ii,!!! 

• BR{Bu T+v) ^ 



AMbo and AM^o [26|, |29| . For the hadronic parameters we follow [29| : 



P^R = -0.71 , P^'^ = -0.9 , Pi 
The remaining parameters used are given in table [T] 
5.2. Lepton sector 

In the Icptonic sector a similar strategy is used: all parameters are evolved to mz and then running 
masses and mixing matrices are used as input for the observables. The implemented formulas are based 
on 



1>LR 



SLL 



= -0.37 , P, 



SLL 



= -0.72 . 



SUSY contributions to the anomalous magnetic moment of the leptons 31 1 



• electric dipole moments (EDMs) of the leptons [32 

• two body decays /i 67, t — > 67 and t 



Ml33| 



• three body decays fi — > ee+e , t — ee+e and r — > fifi^ fJ. 35| 

• Z decays, Z -J> e^/x"^, Z ^ e±T=F and Z ^^r^, [3^ 

5.5. Other constraints 

In addition the EDM of the neutron can be calculated using two different models for the neutron 
where the formulas are based on (sij and we use the same hadronic parameters. Moreover, one can also 
calculate the SUSY contributions to the p-parameter as given in js^]- 



6. Extensions to SLHA 

In this section we describe the SPheno specific extensions to the SUSY Les Houches Accord (SLHA) 
0, Q • We start first with extensions to existing blocks and then discuss new blocks which either control 
the behaviour of SPheno or contain additional model parameters for MSSM extensions. Note, that all 
additional Yukawa couplings have been implemented in complex forms and the corresponding information 
can be passed by using the corresponding blocks starting with IM ^] . 
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6.1. Extensions of existing blocks 

6.1.1. BlockmU'PkR 

In the case that generation mixing is switched on, i.e. the entry 6 contains a non-zero value, then 
independent of this value flavour violation is switched on in the (s)lepton as well as in the (s)quark sector. 

In case of extending the model by a minimal SU{5) as used in fisl] this block gets extended by the 
following entries 

7: 5*0(10) scale where the universal soft SUSY breaking parameters are defined. 
8: extra D-terms due to the breaking of 5*0(10) to SU{5) 
9: A-coupling of the Higgs 24-plet to the 5h 
10: A'-coupling of the Higgs 24-plet to the 5h 

6.1.2. Block MODSEL 

Four switches have been added to flag 3 (particle content) 

2: includes the particle content of a minimal SU{h) model between Mqut and a user chosen 50(10) 
scale, where the SUSY boundary conditions are set. The details of this model are described in 
[l5| . In this case the mass parameters of the right handed neutrinos are stored in the block MNURIN 
(section [6.2.2[) and the corresponding neutrino Yukawa couplings can be stored in the block YNUOIN 
(section 16.2.71) . The data is understood to be defined at the GUT-scale. The additional SU{5) 
parameters as well as the 5*0(10) scale are specified as extensions of the block MINPAR, see section 

3: includes three right-handed (s)neutrinos with a common mass for all three neutrinos. The neutrino 
Yukawa couplings can be specified at the GUT-scale, see section 16.2.71 and the mass of the 
right-neutrinos at their proper scale, see section 15.2.21 

4: includes three right-handed (s)neutrinos which are included at their proper mass scale. The neutrino 
Yukawa couplings can be specified at the GUT-scale, see section 16.2.71 and the masses of the 
right-neutrinos at their proper scale, see section IB. 2. 2 1 

5: includes a Higgs triplet (15-plet) to realize the seesaw II where the formulas of [l3| including the 
corrections presented in [ll| and the 2-loop contributions to the RGEs of the gauge couplings and 
gaugino mass parameters have been implemented. The additional model data are specified in the 
blocks Higgs3IN and YTO, see sections 16.2.11 and I6.2.8l respectivelv. 

10: includes three Higgs 24-plets to realize the seesaw type HI where the complete 2-loop RGEs as 
given m are used. The additional model data are specified in the blocks MWMIN and YB3IN, see 
sections 16.2.31 and 16.2.61 respectively. 

11: includes a Higgs 15-plet to realize the seesaw II where the complete 2-loop RGEs as given in [l^ 
are used. The additional model data are specified in the blocks HiggsSIN and YTIN, see sections 
I6.2.1l and l6.2.9l respectively. 

6.2. New input blocks 

In the output the blocks will be given without the ending IN. 

6.2.1. Block HiggsSIN 

Used to specifiy the model of the seesaw model type II. The data are given in the format 

(2x,i3,2x,lp,el6.8,0p,2x, '# ',a) 

and the entries correspond to 

1: mass of the Higgs triplet (15-plet) 
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Table 2: Default values for fitting R-parity violating parameters if the entries in block NeutrinoBoundsIn are not specified. 
The values are taken from [38| and correspond to the 1 a range but for |(7e3, max P which is 90% CL. 



tan^ Oatm,min 0.8182 
t&V? 6atm,max 1.3256 


tan^ esol,n^in 0.4286 
tan2 Osoi.max 0.4970 




\UeZ^,max? 0.035 


Am^,„^_„ 2.36 • 10-^ cV^ 
^ml,^,ma. 2.54 • 10-3 eV2 


Am^„,_,„,„ 7.46 • 10-5 eV^ 
Am^,™,, 7.83 • 10-5 eV^ 





2: real part of Ai 

3 imaginary part of Ai 
4: real part of A2 

5 imaginary part of A2 

6: if (1) use RGEs for the triplet (15-plet) case 

6.2.2. Blockmmi]i 

In this block one can specify the masses of the right-handed neutrinos within the seesaw I model. The 
masses mm are specified in the FORTRAN format 

(Ix,i3,3x,lp,el6.8,3x, '#' ,a) . 

6.2.3. Block MWMIN 

Here one can specify the mass matrix of the 24-plets Mwij at A/qut for the seesaw type III model 
using the formulas of II3], where the data are given in the FORTRAN format 

(Ix,2i3,3x,lp,el6.8,3x, ,a) 

where the first two integers in the format correspond to i and j and the double precision number to the 
mass parameter. 

6.2.4. Block NeutrinoBoundsIn 

One can use SPheno to obtain R-parity violating parameters consistent with neutrino data. The 
corresponding default values are given in table 16.2.41 This block can be used to modify them. The 
FORTRAN format is 

(Ix , i2 , 3x , Ip , e 16 . 8 , Op , 3x , # , Ix , a) } 

and the entries correspond to 

1: Ato^j„j . . . lower bound on the athmospheric mass difference 
2: Arn^j„j . . . upper bound on the athmospheric mass difference 
3: tan^ Oatm,min ■ ■ ■ lowcr bound on the tan squared of the athmospheric mixing angle 
4: tan^ datm,max ■ ■ ■ Upper bound on the tan squared of the athmospheric mixing angle 
5: ^ra^^i . . . lower bound on the solar mass difference 
6: ^rri^^i . . . upper bound on the solar mass difference 
7: tan^ Osoi,min ■ ■ ■ lower bound on the tan squared of the solar mixing angle 
8: tan^ Osoi.max ■ ■ ■ uppcr bound on the tan squared of the solar mixing angle 
9: C/g3 ,„j„ . . .lower bound on the mixing element Ues squared (reactor angle) 
10: C/g3 ^^j, . . . upper bound on the mixing element Ue3 squared 
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6.2.5. Block SPheno Input 

This block sets the SPheno specific fiags. The FORTRAN format is 

( Ix , i2 , 3x , Ip , el6 . 8 , Op , 3x , # , Ix , a) } 

and the entries correspond to 
1: sets the error level 

2: if 1 the the SPA conventions js^ are used 

3: takes a spectrum which is given by an external program 

4: introduces an extension of the SLHA output: in the case of flavour violation, fiavour ordered states 
are used instead of mass ordered states. 

11: if 1 then the branching ratios of the SUSY and Higgs particles are calculated, if then this calcu- 
lation is omitted. 

12: sets minimum value for a branching ratios, so that it appears in the output 

21: if 1 then the cross sections of SUSY and Higgs particles in e+e" annihilation are calculated, if 
then this calculation is omitted. 

22: sets the center of mass energy Ecms, can be repeated up to 100 times 

23: sets the electron polarisation Pm, can be repeated up to 100 times 

24: sets the positron polarisation Pp, can be repeated up to 100 times 

25: whether to use initial state radation in the calculation of the cross sections 

26: sets minimum value for a cross section, so that it appears in the output 

31: sets the value of Mqut, otherwise Mqjjt is determined by the condition gi — g2 

32: sets strict unification, i.e. gi = 52 ~ 93 

34: sets the relative precision with which the masses are calculated, default is 10~^ 

35: sets the maximal number of iterations in the calculation of the masses, default is 40 

36: whether to write out debug information for the loop calculations 

38: this entry sets the loop order of the RGBs: either 1 or 2, default is 2, i.e. using 2-loop RGBs 

41: sets the width of the Z-boson Tz, default is 2.49 GeV 

42: sets the width of the W-boson Tw, default is 2.06 GeV 

80: if not set the program exists with a non-zero value if a problem has occured 

90: if 1 add R-parity to a high scale spectrum calculated either from mSUGRA, GMSB or AMSB 
boundary conditions 

91: if 1 than bilinear parameters are calculated such that neutrino data are fitted in the experimenatal 
allowed range (the range can be changed using the Block NeutrinoBoundsIn, see section IB. 2. 4p 

92: if 1 gives in case of R-parity violation only the 4x4 MSSM part of the neutrino/neutralino mixing 
matrix N and the correspondingly the 2x2 parts of the charged lepton/chargino mixing matrices 
U and V as well as the block for the stau mixing. This is in particular useful in case one uses the 
program Prospino (40| or older versions of the program Phythia [41| . 
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6.2.6. Block YB3m 

Here one can specify the neutrino Yukawa YiY^ coupling at Mgut for the seesaw type III model using 
the formulas of [13[, where the data are given in the FORTRAN format 

(Ix,2i3,3x,lp,el6.8,3x, '#' ,a) 

where the first two integers in the format correspond to i and j and the double precision number to 
Yukawa coupling. 

6.2.7. Block YNUOIN 

This block specifies the neutrino Yukawa couplings at the GUT scale and the corresponding 
superpotential term is given by = Y^^ijv'^ LjHu. It is assumed that the right-handed neutrinos are in 
the mass eigenbasis. The real parts are specified in the block YNUOIN with the FORTRAN format 

(Ix,i2,lx,i2,3x,lp,el6.8,0p,3x, '#' ,lx,a) 

and the imaginary parts in the block IMYNUOIN with the same FORTRAN input. 

6.2.8. Block YTOIN 

Here one can specify the neutrino Yukawa Y^ coupling at Mqut for the seesaw type II model using 
the formulas of [lO, Eli- The data is given in the FORTRAN format 

(Ix,2i3,3x,lp,el6.8,3x, '#' ,a) 

where the first two integers in the format correspond to i and j and the double precision number to 
Yukawa coupling. 

6.2.9. Block YTIN 

Here one can specify the neutrino Yukawa Y^ coupling at Mqut for the seesaw type II model using 
the formulas of [13[, where the data is given in the FORTRAN format 

(Ix,2i3,3x,lp,el6.8,3x, '#' ,a) 

where the first two integers in the format correspond to i and j and the double precision number to 
Yukawa coupling. 

6.3. New output block 

6.3.1. Block SPhenoLowEnergy 

In this block the calculated values of the low energy observables are given: 

1 BR{b S7) 

2 BR{b sfi+fi-) 

3 BR{b ^ sY^^v.Vi) 

4 BR{B°^ fi+^i-) 

5 BR{B° fi+^i-) 

6 BR{Bu T+v) 

7 BR{Bu ^ T+u)IBR{Bu ^ t+v)sm 

8 A(AfBo) [in ps^i] 

9 A(AfBo) [in ps-i] 

20 SUSY contribution to the anomalous magnetic moment of the electron A(^^)e 

21 SUSY contribution to the anomalous magnetic moment of the muon A(^^)^ 
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22 SUSY contribution to the anomalous magnetic moment of the tau A(^^)t- 

23 electric dipole moment of the electron de 

24 electric dipole moment of the muon 

25 electric dipole moment of the tau dr 

26 BR{n 67) 

27 BR{t e-/) 

28 BR{t fi-f) 

29 BR{^i+ e+e+e") 

30 BRiT+ e+e+e--) 

31 BR{t+ IJ.+li+IJi~) 

39 SUSY contribution to the p-parameter 

40 BR{Z^ e+ix^) 

41 BR{Z^ e+T^) 

42 BR{Z^ ij.+t'^) 

7. Installation and implementing new models 

7.1. Installation 

SPheno can be downloaded from 

http: //projects .hepforge . org/spheno/ 

where the latest tar-ball SPheno3.x.y.tar .gz can found as well as older versions. Unpacking will create 
the directory SPhenoS . x . y where x and y are integers corresponding to the sub- version. This directory 
will contain the following subdirectories: 

• bin: here the executable SPheno will be stored 

• doc: contains the SPheno documentations 

• include: here all the mod- files are stored 

• input: contains input example files 

• lib: here the library libSPheno . a will be stored 

• output: contains the output files corresponding to the examples stored in input 

• src: contains the source code 

The directory SPhenoS. x.y contains a Makefile which can be used to compile SPheno. The default 

compiler is Intels ifort, but by typing make F90=compiler on the console one can use a different compiler 
where compiler has to replaced by the compiler's name. The following compilers have been added NAG 
nagfor, Lahey lf95 and g95. 

It is well known that compilation of the module RGEs . F90 can be time consuming due to the length 
of the 2-loop RGEs for the seesaw models of type II and type III. For this reason they are not compiled 
by default. If the corresponding RGEs should be included then the line 

PreDef = -DGENERATIONMIXING -DONLYDOUBLE 

should be replaced by 
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PreDef = -DGENERATIONMIXING -DONLYDDUBLE -DSEESAWIII 
i.e. add -DSEESAWIII. 

In the case that one want to have quadrupole precision in various parts of the code instead of double 
precision, one has to take out the -DONLYDDUBLE in the Une mentioned above. Note that this can 
substantiahy slow down SPheno. Moreover, not all parts are yet implemented with quadrupole precision. 
The main focus has been on the loop functions as well as on mixing between neutralinos and neutrinos 
in case of R-parity violation. 



1.2. Implementing new models 

New models can easily implemented using the SARAH package 
to put the code generated by SARAH in a new directory within the directory SPhenoS . x . y and run the 
corresponding Makefile. An additional executable will be stored in the directory bin. 



421 . |43|. For this purpose one has 



8. Input and output 

Starting with version SPheno 3.1 there are two main differences with respect to the input and output 

1. SPheno accepts only the SLHA input format as specified and all the output is given in this for- 
mat. In section [5] we have described the extensions to control program specific features as well 
as model extensions. The orginal SPheno input using the files HighScale . in, StandardModel . in 
and Control . in as well as the output in the file SPheno . out have been disabled. Detailed error 
messages and warnings will also be written to the file Messages . out. 

2. One can provide input name and output name as command line options where the first (second) 
name, if present, is interpreted as input (output) filename, e.g. 

SPheno InName OutName 
takes InName for the file containing the input and will write the output to the file OutName. In case 
that the file InName is not found SPheno will look for a file called LesHouches . in as default. The 
default name for the output is SPheno. spc. The length of the names InName and OutName must 
not exceed 60 characters. 



9. Conclusions and comments 

SPheno is constantly developing, in particular in view of implementing additional models and low 
energy observables. In addition it is planed 

• to implement the missing pieces of the SLHA conventions as listed in appendix [Appendix B 

• mixing between Aq and Hq in case of CP phases 

• low energy observables for the case of R-parity violation 

In section [S] several hadronic parameters for the calculation of low energy observables are hard-coded in 
the program. It is planned to construct routines to allow user defined changes in the future. 
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Appendix A. Default SM values 

The following default values will be used if not given in the file LesHouches . in. 

• CKM-matrix, Wolfenstein parameters: A = 0.2265, A = 0.807, p = 0.141, 77 = 0.343 

• gauge sector: l/aem(0) = 137.0359895, mz = 91.187 GeV, Gp = 1.16637 • lO-^GeV-^, 
aY^{mz) = 0.1184 

• lepton masses: = 510.99891 keV, = 105.658 MeV, rrir = 1.7768 GeV 

• quark masses: m„(2 GeV) = 3 MeV, m<j(2 GeV) = 5 MeV, m,.(2 GeV) = 105 MeV, m^imc) = 
1.27 GeV, mi,{mi,) = 4.2 GeV, mt = 171.3 GeV; the top mass is interpreted as on-shell mass 

Appendix B. Unsupported SLHA features 

Here we list the features of the SLHA conventions 0, 0| which are not yet supported: 

• In Block MDDSEL the following entries are currently ignored: 

11: the possiblity to give the parameters at n equidistant values of Q between mz and Qewsb 
21: to give the mass parameters at the pole mass of each individual particle 

• In Block EXTPAR the following entries are currently ignored: 

27: pole mass of the charged Higgs boson 
51: (GVISB only) U{1)y messenger index 
52: (GMSB only) SU{2)l messenger index 
53: (GMSB only) SU{3)c messenger index 



• 


the 


Block 


QEXTPAR 


• 


the 


Block 


RVLAMLLEIN 


• 


the 


Block 


RVLAMLQDIN 


• 


the 


Block 


RVLAMUDDIN 


• 


the 


Block 


RVTLLEIN 


• 


the 


Block 


RVTLQDIN 


• 


the 


Block 


RVTUDDIN 


• 


the 


Block 


RVDIN 


• 


the 


Block 


RVM2LH1IN 



These features will be implemented within the next updates. 

Appendix C. Error messages and warnings, interpretation of the variable kont 

Here we describe how to interpret the values of the variable kont which is used in the error system of 
SPheno. The corresponding warnings and error messages are also given in the file 'Messages. out' if the 
error level is set to the appropriate value. 
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Appendix C.l. Module Mathematics 
-1: step size gets to small in routine ODEint 

-2: maximal value > 10^'' ODEint 

-3: too many steps are required in routine ODEint 

-4: boundary conditions cannot fuUfilled in routine ODEintB 

-5: maximal value > 10^^ ODEintB 

-6: step size gets too small in routine ODEintB 

-7: too many steps are required in routine ODEintB 

-8: boundary conditions cannot fullfilled in routine ODEintC 

-9: maximal value > lO''^'' GDEintC 
-10: step size gets too small in routine ODEintC 
-11: too many steps are required in routine ODEintC 
-12: step size gets too small in routine rkqs 

-13: the size of the arrays do not match in routine ComplexEigenSystems 

-14: potential numerical problems in routine ComplexEigenSystems 

-15: the size of the arrays do not match in routine RealEigenSystemse 

-16: potential numerical problems in routine RealEigenSystems 

-17: the size of the arrays do not match in routine tqli 

-18: too many iterations in routine tqli 

-19: too high accuracy required in routine Dgauss 

-20: too high accuracy required in routine Dgauss Int 

-21: precision problem in routine Kappa 

-22: step size gets too small in routine IntRomb 

-23: too many steps are required in routine IntRomb 

-24: singular matrix in routine GaussJ 

-25: inversion failed in routine InvMatS 

-26: singular matrix in routine GaussJ 

Appendix C.2. Module StandardModel 

-101: routine CalculateRunningMasses: Qiow > rnb{mb) 

-102: routine CalculateRunningMasses: Max{Qiow,'rni,{mh) > Qmax) 
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Appendix C.3. Module SusyMasses 

-201: negative mass squared in routine ChargedScalarMassEpslnt 

-202: negative mass squared in routine ChargedScalarMassEps3nt 

-204: \Yr\'^ < in routine CharginoMassS 

-205: \Yt\^ < in routine CharginoMassS 

-206: negative mass squared in routine PseudoScalarMassEpslnt 

-207: negative mass squared in routine PseudoScalarMassEpsSnt 

-208: negative mass squared in routine PseudoScalarMassMSSMnt 

-210: negative mass squared in routine ScalarMassEpslnt 

-211: negative mass squared in routine ScalarMassEps3nt 

-212: negative mass squared in routine ScalarMassMSSMef f 

-213: negative mass squared in routine ScalarMassMSSMnt 

-215: m|o < in routine ScalarMassMSSMef f 

-216: m^o < in routine ScalarMassMSSMef f 

-217: m|+ < in routine ScalarMassMSSMef f 

-220: negative mass squared in routine Sf ermionMasslEpsl 

-221: negative mass squared in routine Sf ermionMasslEps3 

-222: negative mass squared in routine Sf ermionMasslMSSM 

-223: negative mass squared in routine Sf ermionMassSMSSM 

-224: negative mass squared in routine SquarkMassSEps 

-225: m| < in routine TreeMassesEpsl 

-226: m| < in routine TreeMassesMSSM 

-227: m\o < in routine TreeMassesMSSM 

-228: < in routine TreeMassesMSSM 

-229: m? < in routine TreeMassesMSSM2 

-230: m\o < in routine TreeMassesMSSM2 

-231: mjj+ < in routine TreeMassesMSSM2 

-232: m? < in routine TreeMassesMSSM3 
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Appendix C.4- Module InputOutput 

-302: routine LesHouches_Input: unknown entry for Block MODSEL 

-303: routine LesHouchesJEnput: model must be specified before parameters 

-304: routine LesHouchesJEnput: unknown entry for Block MINPAR 

-305: routine LesHouchesJEnput: model has not been specified completly 

-306: routine LesHouchesJEnput: a serious error has been part of the input 

-307: routine LesHouchesJEnput: Higgs sector has not been fully specified 

-308: routine ReadMatrixC: indices exceed the given boundaries 

-309: routine ReadMatrixR: indices exceed the given boundaries 

-310: routine ReadVectorC: index exceeds the given boundaries 

-311: routine ReadVectorR: index exceeds the given boundaries 

-312: routine ReadMatrixC: indices exceed the given boundaries 

Appendix C.5. Module SugraRuns 

-401: routine BoundaryEW: negative scalar mass squared as input 
-402: routine BoundaryEW: m%{mz) < 
-403: routine BoundaryEW: sin^ 0-^ < 
-404: routine BoundaryEW: m'^ < 

-405: routine BoundaryEW: either m/oij/m; < 0.1 or mi^B./mi > 10 
-406: routine BoundaryEW: cither md^R/mu < 0.1 or mdoR/md > 10 
-407: routine BoundaryEW: cither mu^u/md < 0.1 or ?Ti„j-,/{/m„ > 10 
-408: routine RunRGE; entering non-perturbative regime 
-409: routine RunRGE: nor gi ^ 52 at Mqut neither any other unification 
-410: routine RunRGE: entering non-perturbative regime at Mgut 
-411: routine RunRGE: entering non-perturbative regime at Mh^ 
-412: routine Sugra: run did not converge 
-413: routine Calculate_Gi_Yi: m%{mz) < 

-414: routine Calculate_Gi_Yi: too many iterations to calculate mb{mb) in the MS scheme 
-415: routine Sugra: < at niz 
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Appendix C.6. Module LoopMasses 

-501 negative mass squared in routine SleptonMass_lL 

-502 p^ iteration did not converge in routine SleptonMass_lL 

-503 negative mass squared in routine SneutrinoMass_lL 

-504 p"^ iteration did not converge in routine SneutrinoMass_lL 

-505 negative mass squared in routine SquarkMass_lL 

-506 p^ iteration did not converge in routine SquarkMass_lL 

-507 m^o < in routine LoopMassesMSSM 

-508 m\o < in routine LoopMassesMSSM 

-509 < in routine LoopMassesMSSM 

-510 |/u|^ > 10^° in routine LoopMassesMSSM 

-511 l/Up < in routine LoopMassesMSSM 

-512 m%{mz)'^ < in routine LoopMassesMSSM 

-513 m^o < in routine LoopMassesMSSM_2 

-514 m\o < in routine LoopMassesMSSM_2 

-515 < in routine LoopMassesMSSM_2 

-516 |;u|2 > 10^° in routine LoopMassesMSSM^ 

-517 l/ip < in routine LoopMassesMSSM^ 

-518 m|(mz)^ < in routine LoopMassesMSSM_2 

-519 m^o < in routine LoopMassesMSSM_3 

-520 m\o < in routine LoopMassesMSSM^ 

-521 < in routine LoopMassesMSSM^ 

-522 l/ip > 10^° in routine LoopMassesMSSM^ 

-523 l/ip < in routine LoopMassesMSSM^ 

-524 m|(mz)^ < in routine LoopMassesMSSM J 

Appendix C.7. Module TwoLoopHigg.sMa.ss 
-601: routine PiPseudoScalar2: m? < 

-602: routine PiPseudoScalar2: m? < 

b 

-603: routine PiPseudoScalar2: m| < 
-604: routine PiScalar2: m? < 
-605: routine PiScalar2: m? < 



-606: routine PiScalar2: m? < 
-607: routine Two_Loop_Tadpoles: m| < 
-608: routine TwoXoop_Tadpoles: m? < 
-609: routine Two_Loop_Tadpoles: m? < 
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Appendix C.8. Module MathematicsQP 
-1001: the size of the arrays do not match in routine ComplexEigenSystems_DP 

-1002: potential numerical problems in routine ComplexEigenSystems_DP 

-1003: the size of the arrays do not match in routine ComplexEigenSystems_QP 

-1004: potential numerical problems in routine ComplexEigenSystemsjQP 

-1005: the size of the arrays do not match in routine RealEigenSystems_DP 

-1006: potential numerical problems in routine RealEigenSystemsJDP 

-1007: the size of the arrays do not match in routine RealEigenSystemsJQP 

-1008: the size of the arrays do not match in routine Tqli_QP 

-1009: too many iterations in routine Tqli_QP 

-1010: too many iterations in routine Tql2_QP 

Appendix D. Loop corrections 

Here we list the improvements which have been implemented in SPheno with respect to ref. 
• in the 1-loop corrections to the gluino mass we use for the gluon contribution 



A(E^) = {B,{p^my,0) - 2i?l(p^m|^,0)) (D.l) 

where to^.t is the tree level gluino mass and which reduces for p^ — m? to the formula 



9l ' ' 



of ref. 0. 



• In addition flavour violation has been taking into account and the corresponding formulas can be 
found in 
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